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Abstract 



A new mechanism for understanding small neutrino masses using only simple 
^ ■ new physics at the TeV scale is proposed. As an application, it is shown how it 

If^ . can naturally lead to the mass hierarchy of the so called bimaximal mixing in 

, the case of three active neutrinos, or the (3+1) scenarios for sterile neutrinos, 

I using only the SU{2)l quantum numbers of the particles. 
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I. INTRODUCTION 



The most commonly used way to understand the smallness of neutrino mass compared 
to that of the charged fermions is the seesaw mechanism [Q], according to which one has, 

_^ : -M. ^ -MlM],'Mn. (1) 

For choices of the Dirac mass matrix Mr, dictated by simple unified models, one expects 
the typical scale of the seesaw masses to be around 10^° GeV or higher for neutrino masses 
required to fit atmospheric and solar data 0,^. In this class of models, neutrino masses and 
mixings are dictated by physics at superheavy scales. An important question to ask (and 
has often been discussed in literature) therefore is whether the neutrino masses can arise 
from new physics nearby the weak scale, e.g., in the TeV range so that one may hope to 
probe them in the collider or other indirect experiments. 

There exist several TeV scale scenarios for neutrino masses in the literature: (i) radiative 
scenarios that use extensions of the standard model by extra charged Higgs bosons P] (ii) 
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scenarios using bulk neutrinos or bulk scalars P,[l0| and (iii) a more recent one |TT 
that uses the seesaw mechanism but with naturally suppressed Dirac masses M£,. 



In this letter, we propose a new mechanism that uses higher dimensional operators |T2 
that lead to a naturally suppressed neutrino masses using only simple new physics at the 
TeV scale. The higher dimensional operator could arise in a number of ways from TeV scale 
physics depending on the nature of extra symmetries and the particle content of the theory. 
We will illustrate our mechanism using the effective low energy theory at the weak scale. 

For simplicity, we will assume the low energy theory to be the standard model plus an 
extra Higgs doublet and a Higgs singlet field. We will assume also the theory to have softly 
broken global symmetries such as L = Lg + + L,- and other linear combination of these 
quantum numbers depending on the kind of neutrino mixings that one wants to obtain. The 
general principle we will adopt is that global symmetries will only be broken by soft terms 
in the potential and any other interaction with dimension d > 4 will respect the symmetry. 

As an application of our mechanism, we apply it to understand the lightness of the sterile 
neutrino, that is required to provide a simultaneous fit to solar, atmospheric and the LSND 



data |]T3[. We find that in simple models that implement our strategy, the sterile neutrino, 
even though is a standard model gauge singlet, may get a small mass comparable to that of 
familiar neutrinos due to the higher dimensional operators, leading naturally, for instance, to 



a (3+1) scenario which has recently been discussed [1^,15|. We must however mention that 



the (3+1) scheme is generically constrained such that it may work only for certain domain 
of values of the mass difference ^m\g^Q [|14|,|15|, and there are also recent papers [|1^ noting 
that it may be lest favoured as the explanation of all neutrino anomalies, though it has not 
been ruled out, and so, it remains of some interest. 



II. THE NEW MECHANISM 

Before presenting our proposal, let us review the higher dimensional mechanism suggested 



by Weinberg [12]. According to this mechanism, one writes down the nonrenormalizable 



operator involving the standard model fields that can lead to neutrino masses, 

fab ( jT n-i 



^ — [LiaC LjbHkHi£ikejij , (2) 

where i,j,k,l are SU{2)l indices and A is the scale of new physics. This operator breaks 
explicitly the lepton number. When H develops a vacuum expectation value (vev), neutrinos 
acquire masses, 

= /i^. (3) 

From Eq. (^), it is clear that this formula is similar to the seesaw formula in Eq. (1) with 
Dirac mass replaced by the Higgs vev that determines the weak scale and to get neutrino 
masses in eV range or less implies the scale of new physics A > 10^^ GeV or so. We should 
mention that a recent study [|l7l has classified several other non renormalizable operators of 



dimension smaller than 11 that may give rise to neutrino masses. However those equivalent 
to the ones we will use here were not considered there. 
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Let us now add to the SM representation content a singlet scalar field x which carries 
lepton number L = — 1 and a second scalar doublet H' with L = 0. We then impose the 
discrete Z2 symmetry on the model under which H' and x odd while all other fields 
are even. Then the higher dimensional operator given above is forbidden by lepton number 
conservation. The lowest dimensional operator that is invariant under L is given by: 

^-Jf,{LHxf^{^,{LH'x)\ (4) 

where M is now the new scale of physics. Once H, H' and x develop vev's we obtain the 
following expression to neutrino masses: 

^^-^^{fW' + f{H'f). (5) 

If the vev of (x) is suppressed, then M can be in the TeV range. We now show that this 
happens quite naturally in our model. 

To study the vevs of the Higgs fields, we write down the most general Higgs potential 
consistent with the symmetries we have imposed on the model, i.e., the gauge symmetry and 
global U{1)l. However to generate a natural small vacuum to x we consider in the potential 
terms that break explicitly the lepton number, 

V{H, H', x) = iAH^H + ^{H^Hf + iJi]j,H'^H' + ^{H'^H'f + X^H^HH'^H' 

+^X^X + y(x^x)' + KH^Hx^x + KH'^H'x^x - MH^^H'x + h.c. (6) 

First note that the M term in the potential breaks the global L-symmetry softly. For 
n'jj Ht < 0, the potential leads to both H and H' fields having nonzero vev's that break the 
electroweak symmetry. Minimization of the potential gives for the parameterization. 





the mixing angle: 



(7) 



If one assumes the dn and en to be odd under the Z2 symmetry, then one can choose 
{H') <^ (H) as in the case of large tan j3 supersymmetric models. For /x^ > and A2 > 0, 
the M term in the potential then induces a vev for the % field given by: 

^ M{H){H') ^ Mvuv, 

M2 - /x2+A^(x^) + ^;2cos/?2(A4 + A^tan2/?)- 

If we assume ~ and ~ M ~ 20 TeV, we get (x) ^0.3 GeV, without any fine 
tuning of parameters. Note that this leads to breaking of lepton number. There is however 
no zero mass particle since the M term in the potential breaks lepton number explicitly. In 
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fact the masses of the real fields in x become of order of the the scale M. It is worth noticing 
that this is nothing but a singlet realization of the so called Type II see-saw mechanism 



Indeed in all our analysis above one can use a electroweak triplet carrying the same lepton 
number instead of the singlet Xi with the richness on the phenomenology of such a field. For 
what we want to discuss hereafter, whether x is a singlet or a triplet will not be relevant. 

Using Eq. and / ~ 0.1, we get ~ 0.04 eV. It is interesting that the muon neutrino 
mass required to explain the atmospheric neutrino deficit emerges without any fine tuning 
of parameters. This is the new mechanism which we want to exploit in the rest of the paper 
to understand neutrino mass patterns. This mechanism has also the potential to be useful in 
understanding small neutrino masses in low string scale theories, where the highest available 
scale is in the multi TeV range. Moreover, it may even generate a naturally small sterile 
neutrino mass. We should notice that the present mechanism can also be seen as a Type II 
see-saw realization of the dimension nine operator L^H^. Notice that such class of operators 
was explicitly excluded in the general analysis of Ref. |T7| . 

Let us conclude this section describing an underlying theory that could lead to the higher 
dimensional operators as part of the low energy effective theory. One possibility is to have at 
the 10 TeV scale three new singlet fields of mass of order 10 TeV: such that have 

L = ±1 respectively whereas has L = and is odd under Z2- The invariant potential 
involving these fields is then given by: 

C = LHN2 + M1N1N2 + NsxNi + M3N3N3 + h.c. (10) 

Ml 3 are assumed to be of order 10 TeV. This theory leads in the low energy limit to the 
effective theory described in (H) (see Fig. 1). We now move into considering how the desired 
pattern for the neutrino masses can be obtained. 



III. BIMAXIMAL MIXING PATTERN 



In this section, we apply the mechanism of the previous section to derive the bimaximal 
neutrino mixing pattern. For this purpose, we recall that in the previous section, we used 
the softly broken global symmetry L. Suppose, in addition, we have a softly broken global 
V = Lf. — — Lt-. Then, the allowed higher dimensional terms in the Eq. (|^) would have 
generation labels in them and would lead to a mass matrix of the form: 



/ ce se 
ce 
V^e 



(11) 



where cg [sg) represents the function cos^ (sin 6*). It is well known that this leads to the 
bimaximal neutrino mixing pattern (see for instance Refs. p!9|-p^): 



U 



BM 



C0 

\ V2 



1 

V2 



\ 



(12) 



As it is evident from ([TTD , at this stage, there is no mass splitting between ui and z/2 
that can lead to oscillations. There is however mass splitting between ui 2 and z/3 with 
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Am^ = m^, which give rise to oscillations with the amplitude sin^2^. As we saw before, 
for the parameters chosen in the previous section, this is precisely in the range required to 
solve the atmospheric neutrino anomaly, since then Am^j^ = 1.6 x 10~^ eV^. The 

interesting point to reemphasize in this connection is that the neutrino mass spectrum is an 
inverted one. 

In order to generate the splitting between the ui and 1^2, we need to introduce soft 
breaking of L' = — L^^ — Lr symmetry so that an entry such as ruee or m^^ can be 
generated. For solar neutrino anomaly solution via the large angle MSW mechanism, one 
would require to be of order 10~^ eV and for resolution via the vacuum oscillation 
solution, one would require m^e ~ 10~^ eV. For this purpose it is convenient to introduce 
a field 5 with L = but L' = —1. The vev of this field can also be induced by the 
mechanism developed in the last section, but this time we choose the parameters such that 
{S)/M ~ 10~^. This allows us to introduce a higher dimensional operator of the type: 

{L,Hfx'6'/M'. (13) 

Once all vevs are substituted, one gets rriee ~ 4 x 10~^ eV. This favors the large angle MSW 
solution to the solar neutrino problem. In fact, the solar mass parameter is then given as 
Ami ~ 2m2((5)/M)2 ^ S-IQ-^ eV^ A smaller vacuum, (6) ~ 1 GeV, however, will produce 
rripp ~ 3 X 10~^ eV for vacuum solution. 



IV. THE (3+1) SCHEME 

The above neutrino mass scheme does not explain the LSND results. In order to include 
it into our model, the simplest possibility in our framework appears to be the so called 



(3+1) scheme, that may account for the explanation of all neutrino anomalies [p!4|-[T6[|. For 
this purpose, we include a left-handed singlet neutrino in the model, Ug, that carries its 
own lepton number Lg, and two odd scalar singlets, besides X; which we call 61 and 62. 
In order to get the final mass pattern, extra softly broken global symmetries are needed. 
Following our principle, we take those symmetries to be linear combinations of the four lepton 
numbers: Le^^,T,s- Let us consider, for instance, three of such global symmetries: U{1)li', 
[7(1)^2 and [7(1)^3, where i^i,2,3 represents our specific choice for the linear combinations 
of lepton numbers. We have taken such symmetries among the total active lepton number, 
L = Lf. + Lf^ + Lt] the sterile lepton number, L^; the total lepton number, Lj- = L + L^; 
and the combinations: V = L^. — — Lj. and L" = Lg — L. We always pick up L2 = V 
to maintain as much as possible our previous results. The quantum numbers associated to 
the chosen symmetries and carried by the scalar fields are depicted in Table I, where three 
different models are considered. All these models have the same output, they provide the 
following higher dimensional mass terms (at the lower order) that conserve such symmetries: 
i) Active flavour mass terms: 

J I] ' ^^^^ 

where the Yukawa couplings fij] for i = efir, are only constrained by the L' symmetry. 
These terms give rise to the same mass matrix as in Eq. (|ll|), and same order of parameters, 
that now we assume. 
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ii) Active-sterile couplings: 



IP + ^^^^ + ^"^^ M3 • ^^^^ 

In this case the Yukawa couplings are all non zero, and if (^1,2) are both the same order then 
the mass couplings, fii ~ (H) (x)'^ (S) / , are comparable among each other. For simplicity 
one may take the light hierarchy ge : : gr ^ ^ '■ ce : sg. This choice is suggested by the fact 
that they are governed by L' as well as the active mass terms. Some amount of tuning may 
be needed in our example, but small deviations of our choice will not affect our conclusions. 
However, this particular choice will allow us to make more explicit calculations later on. 
iii) Finally, the sterile mass term: 

That give rise to the sterile mass, rris ~ (x)"^ {Si) {62) / . Now, if we take (x) ~ -3 GeV, as 
before with M ~ 20 TeV and (6) ~ 400 GeV, we get for /i ~ 1, the sterile mass ~ 2 eV, 
which is adequate to explain the mass difference required by the LSND data. 

After introducing expectation values we get, in the basis {ua, ^s), the mass matrix, 

\ e ' rUs ms I 

where ei = Hi/nis ^ {gi/h){H)/ (S). Thus, by taking g ~ C(l), this give us the hierarchy on 
masses, m <^ Hi < rris that will reproduce the features of the (3-1-1) scheme. 

Let us now analyze how the required parameters appear from Eq. (jTZ]). First thing to 
notice is that at the lower order in e, the unitary mixing matrix has the form: 



Therefore, as an immediate conclusion we have that t/gs = 0, as desired by CHOOZ. Also, 
LSND arises with siv? Olsnd = 4|eee^p, that constrains the active to sterile couplings to 
values well in our preferred range of parameters (O(10~^)). 

Next, while the sterile mass remains almost unperturbed, all other mass terms in the 
active block get a see-saw type correction 

^ A<W -e e V . (19) 

The effect of this correction is to break the degeneracy of ui and 1/2 induced by L', since 
in general it introduces small diagonal terms (among others). Thus, a nice feature of this 
model is that we do not need to invoke any other mechanism to generate those corrections. 
They come naturally due to the breaking of the other symmetries. 

Let us now consider as an specific example the above suggested hierarchy of the g Yukawa 
couplings. We then take e ~ (|e|/v^)(l, cg, sg). After rotating the active sector by Ubm, 



A^(°^ becomes Diag{m, —m,0), where m = y Am'^tm ~ 0-04 eV, as before. The same 
rotation projects e along the direction (1,0,0). Hence, z/3 remain exactly massless, and 
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the 1-4 sector remains mixed with a mixing angle tan/3 ~ |e|. Therefore, after a complete 
diagonalization, one gets the exact mixing matrix: 



/ 




-1 


V2 


V2 






C0 




V2 








se 






V2 


v 








^\ 
J 



(20) 



The mass spectrum we then get is, 



mi = m — \e\^ms , 
m2 = —m , 
mg = , 

m4 = m,(l + |ef ) . (21) 

Let us see how the solar neutrino mass arise. Looking at the spectrum (^Tj), we require 
Am2i = rn^ — m\ = Am^,, that renders, 

|ef ^ 7^ ■ 22 

This may imply a quite large cancellation on the parameters of the theory, and is a man- 
ifestation of the constrained solution we are discussing. The required fine tunning would, 
however, be not too large. For large mixing angle solutions one needs that |e|^ be adjusted 
at the level of 10% at most, such that the difference 2m — |e|^ be of the order of 10~^. 

By taking m^ ~ 1 — 2 eV as the LSND scale (just as it came from our previous general 
analysis), for the large mixing MSW solution of solar neutrino problem, one gets: |e|^ ~ 
0.04 — 0.08. Therefore, our prediction for the LSND oscillation probability is 

sin^ 2eLSND = cosl ~ (-8 - 3-2) x 10"^ (23) 
where the extreme right hand side has been evaluated assuming maximal mixing in the 



active sector. We should comment that, as already observed in the literature |T3JTH], this 
also pushes Bugey and CDHS to the limit. With our set of parameters we get: sin^ 26'Bugey ~ 
0.04 — 0.08 and sin^ 2^cdhs ~ 0.08 — 0.16. On the other hand, this could also be interpreted 
as a chance to have accessible signals of Ug^^j^ disappearance in future experiments. Given our 
(inverted) hierarchy, other mixing parameters are given as 

sin^ 2eatm = sin^ 20 ■ (l - ^\e]^^ ^ 0.96 - 0.98 , 
sin^ 2^0 = (1 - lel^) ^ 0.92- 0.96 , (24) 

which means that solar and atmospheric are both explained by large mixing angle solutions 
in our model. 
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V. CONCLUDING REMARKS 



We have presented a new mechanism that may generate small neutrino masses at tree 
level which does not need the existence of large scales for new physics but rather make 
use of higher dimensional operators generated at the TeV scale. In the realization of our 
mechanism the smallness of the neutrino masses are due to the large mass suppression on 
dimension seven operators which are involved in the generation of masses and due to the 
presence of a relatively small vacuum of a scalar particle which (softly) breaks down lepton 
number. Such a vacuum arises from a singlet scalar realization of the so called Type II 
see-saw mechanism. As a direct application we have explored the generation of textures 
that may explain the neutrino anomalies. Given the flavour symmetry L' : ^ — 
we generate the bimaximal neutrino mixing patter, which we further extended by adding 
extra symmetries and a sterile neutrino which also obtains a eV mass, giving rise to the 
(3+1) scenario that may account for including LSND in the scheme. Our mechanism may 
be useful in theories where the fundamental (string) scale is in the TeV range, since it relays 
only on brane physics for the generation of neutrino masses. 
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a Pesquisa do Estado de Sao Paulo (FAPESP). 
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TABLES 





Model I 


Model II 


Model III 


Field 


Lj' 


L' 


L" 


Lj' 


L' 




L 


L' 


Ls 


X 


-1 







1 


-1 










-1 








Si 





-1 




-2 





-1 




-1 


1 


-1 


-1 


S2 





1 




-2 





1 




-1 


1 


-1 


1 



TABLE I. Assignments of charges for the scalar fields in three different models with the flavour 

symmetries as shown, which give rise to the (3+1) scenario. Here Lt = Ls + L; L = Lg + + L-r] 
L' = Le — L^, — L-ri and L" = Lg — L. Standard Higgs doublets are in all cases chargeless under 
this symmetries. 
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FIG. 1. Tree-level realization of the effective operator in Eq. (0). 
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